Partonic Effects on Pion Interferometry at the Relativistic Heavy Ion Collider 
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Using a multiphase transport (AMPT) model that includes both initial partonic and final hadronic 
interactions, we study the pion interferometry at the Relativistic Heavy Ion Collider. We find that 
the two-pion correlation function is sensitive to the magnitude of the parton scattering cross section, 
which controls the parton density at which the transition from the partonic to hadronic matter 
occurs. Also, the emission source of pions is non-Gaussian, leading to source radii that can be 
more than twice larger than the radius parameters extracted from a Gaussian fit to the correlation 
function. 
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Particle interferometry based on the Hanbury-Brown 
Twiss (HBT) effect has long been used to measure the 
size of an emission source jij. In heavy ion collisions, it 
has been suggested that the HBT study can provide in- 
formation not only on the spatial extent of the emission 
source but also on its expansion velocity and emission du- 
ration |^ H In particular, the long emission time 
as a result of the phase transition from the quark-gluon 
plasma to hadronic matter in relativistic heavy ion colli- 
sions is expected to lead to an emission source which has 
a much larger radius in the direction of the total trans- 
verse momentum of detected two particles (i2 ou t) than 
that perpendicular to both this direction and the beam 
direction (i? S idc) [§, ||, @] • Since the quark-gluon plasma is 
expected to be formed in heavy ion collisions at the Rel- 
ativistic Heavy Ion Collider (RHIC), it is thus surprising 
to find that the extracted ratio i? ou t/^sidc from a Gaus- 
sian fit to the measured two-pion correlation function in 
Au+Au collisions is close to one [|[ ||, Also, the 
extracted radius parameters are small compared to the- 
oretical predictions based on the hydrodynamical model 
H . These experimental results, especially the small value 
of i?out/-Rsido, have been attributed to strong space-time 
and momentum correlations in the emission source [ fTT| . 

In this Letter, a multiphase transport model (AMPT) 
p2| , |l3| , |i4[ , that includes both initial partonic and fi- 
nal hadronic interactions, is used to study the pion in- 
terferometry in central Au+Au collisions at RHIC at 
V« = 130 A GeV. The AMPT model is a hybrid model 
that uses the minijet partons from hard processes and the 
strings from soft processes in the HIJING model jL{| as 
the initial conditions for modeling the collision dynam- 
ics. The time evolution of partons is then modeled by the 
ZPC |L6| parton cascade model. At present, this model 
includes only parton-parton elastic scatterings with an 
in-medium cross section given by: 
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where the strong coupling constant a s is taken to be 0.47, 
and s and t are the Mandelstam variables. The effective 
screening mass fi depends on the temperature and density 



of the partonic matter but is taken as a parameter for fix- 
ing the magnitude and angular distribution of the parton 
scattering cross section. After these minijet partons stop 
interacting, they are combined with their parent strings 
to fragment to hadrons using the Lund string fragmenta- 
tion model as implemented in the PYTHIA routine |l7| . 
The final-state hadronic scatterings are then modeled by 
the ART model fl§. 

The default AMPT model has been quite successful in 
describing the measured rapidity distributions of charge 
particles, particle to antiparticle ratios, and the spec- 
tra of low transverse momentum pions and kaons |l3f| 
in heavy ion collisions at the Super Proton Synchrotron 
(SPS) and RHIC. Since the initial energy density in 
Au+Au collisions at RHIC is expected to be much larger 
than the critical energy density at which the hadronic 
matter to quark-gluon plasma transition would occur 
@ |o), the AMPT model has been extended to con- 
vert the initial excited strings into partons |L4|. In this 
string melting scenario, hadrons, that would have been 
produced from string fragmentation, are converted in- 
stead to valence quarks and/or antiquarks. Interactions 
among these partons are again described by the ZPC par- 
ton cascade model. The transition of the partonic matter 
to the hadronic matter is, however, achieved using a sim- 
ple coalescence model, which combines two nearest par- 
tons into mesons and three nearest partons into baryons 
or anti-baryons. Using parton scattering cross sections 
of 6-10 mb, the extended AMPT model is able to repro- 
duce both the centrality and transverse momentum (be- 
low 2 GeV/c) dependence of the elliptic flow measured 
in Au+Au collisions at *Js = 130A GeV at RHIC 

From the AMPT model, the source of emitted parti- 
cles is obtained from their space-time coordinate x and 
momentum p at freezeout, i.e., at their last interactions. 
Denoting the single-particle emission function for pions 
by S(x,p), the HBT correlation function for two identi- 
cal pions in the absence of final-state interactions, such 
as the Coulomb interaction, is then given by pL E2| : 
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where pi and P2 represent the momenta of the two pions, 
respectively, K = (pi + pa)/2, and Q — (pi — p 2 , E\ — 
E-z). Expecting that the emission function is sufficiently 
smooth in the momentum space, one can evaluate the 
correlation function by using pi and p 2 for K in the nu- 
merator of Eq.(||). Furthermore, the three-dimensional 
correlation function in Q is usually shown as a function 
of the invariant relative momentum (Qi nv = \J —Q 2 ) or 
as a function of the projection of Q in the "out-side-long" 
(osl) system J|, Q, defined by the beam direction (Qi on g), 
the direction along the total momentum of the two par- 
ticles in the transverse plane (Qout), and the direction 
orthogonal to the above two directions (Q s idc). 
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FIG. 1: Correlation functions for midrapidity charged pi- 
ons with 125 < pt < 225 MeV/c. Theoretical results with 
Coulomb corrections are shown for the default AMPT model 
and for the extended AMPT model with string melting and 
various values for a v . Coulomb-uncorrected correlation func- 



tions from the STAR collaboration 
circles. 



are shown by filled 



Using the emission function obtained from the AMPT 
model for central (b = fm) Au+Au collisions at y/s 
1 30.1 GeV, we have evaluated the correlation function 
C*2(Q,K) in the longitudinally comoving system frame 
using the program Correlation After Burner |2^] . In Fig. 
|l|, we show the one-dimensional projections of the calcu- 
lated correlation function including final-state Coulomb 
interactions for midrapidity (—0.5 < y < 0.5) charged pi- 
ons with transverse momentum 125 < pt < 225 MeV/c. 
Also shown are the measured ir~ correlation functions 
from central collisions by the STAR collaboration with- 
out removing the effect due to the Coulomb interaction 
||. In evaluating the one-dimensional projections of the 
correlation function onto one of the Q Q ut , Qsidc and Qi on g 



axes, we have integrated the other two Q components 
over the range — 35 MeV/c. The dash-dotted curves 
in Fig. [l] are results from the default AMPT model 
(no string melting) with a parton scattering cross sec- 
tion of cr p = 3 mb, while other curves are those from 
the extended AMPT model with string melting but dif- 
ferent values for a v . It is seen that with string melt- 
ing both the width of the Q mv correlation function and 
its height decrease with increasing a p . The decreasing 
width can be more clearly seen from the calculated cor- 
relation functions without including final-state Coulomb 
interactions [£4j , as the value of the correlation function 
in this case is exactly two at Q mv — 0. Since the integra- 
tion of the other two Q components over a fixed range 
(0 — 35 MeV/c) gives a smaller value when the correla- 
tion function becomes narrower in Q, the height of the 
one-dimensional projection of the correlation function de- 
creases with increasing a p . To reproduce the measured 
one-dimensional correlation functions by the STAR col- 
laboration, we find that a parton scattering cross section 
of about 10 mb is required in the extended AMPT model 
with string melting. Since in our model partons are con- 
verted to hadrons after they make their last scattering, 
the parton scattering cross section thus controls the den- 
sity at which partons are converted to hadrons. A larger 
a p leads to a lower density and larger volume for the 
parton to hadron phase transition. We expect that the 
two-pion correlation data could also be reproduced by 
using a smaller parton scattering cross section but delay- 
ing the hadronization until the critical density is similar 
to the one given by the parton freezeout condition using 
a scattering cross section of 10 mb in the parton cascade 



model |24 



The size of the emission source can be determined from 
the emission function via the curvature of the correlation 
function at Q = 0: 
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= D Xi . Xj {K)-D x ^ jt {K)-D^ X] {K)+D p ^ Pjt {K), (3) 

where Xi(i = 1 — 3) denotes the projections of the particle 
position at freezeout in the osl system, i.e., x ou t, x s id e and 
a^iong, and — K/ifo with Kq being the average energy 
of the two particles. In the second line of Eq.(|J), we have 
the variance D x y — (x ■ y) — (x)(y), with (x) denoting 
the average value of x. 

Since only the correlation function is measured in ex- 
periments, the size of emission source is usually estimated 
by fitting the measured correlation function C2(Q,K), 
that has been corrected for effects due to final-state 
Coulomb interactions, with a four-parameter Gaussian 
function: 



C 2 (Q, K) = 1 + Aexp - Rl(K)Ql 



(4) 
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If the emission source is Gaussian in space-time, then for 
central heavy ion collisions considered here the radius 
parameters obtained from the above Gaussian fit to the 
correlation function would be the same as those deter- 
mined directly from the emission function of the source 
via Eq.(||). However, because of space-time correlations 
in the emission function, such as those induced by the 
collective flow, the radius parameters in general do not 
give the true source size. 
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FIG. 2: Source radii from the emission function (solid curves) 
as well as the fitted radius and A parameters from the Gaus- 
sian fit to the correlation function (curves with squares) for 
midrapidity pions as functions of pion transverse mass ttit- 
The inset shows the corresponding ratio R ou t /-Rsidc • 

In Fig. ||, we show by solid curves the transverse mass 
(wt) dependence of the source radii i? ou t (upper- left 
panel), -R s idc (upper-right panel), and i?i on g (lower-left 
panel) determined from the emission function for midra- 
pidity charged pions given by the AMPT model with 
string melting and parton cross section of a p = 10 mb 
for central Au+Au collisions at y/s — 13QA GeV. The 
radius parameters determined from the Gaussian fit to 
the three-dimensional correlation function in Q without 
including final-state Coulomb interactions are shown by 
curves with squares. For the wit values shown here, these 
radius parameters are about a factor of 2 to 3 smaller 
than the source radii obtained directly from the emission 
function. The emission source from the AMPT model 
thus deviates appreciably from a Gaussian one. The ra- 
dius parameters from the Gaussian fit to the correlation 
function have similar values as the experimental ones 
from a Gaussian fit to the measured correlation function 
after correcting for the final-state Coulomb interactions 
||. The parameter A (scaled up by a factor of 10) from 
the Gaussian fit to the correlation function is shown in 



the lower-right panel by the curve with squares. It has 
a value of about 0.5 at low tot but increases to about 
1 at large tot- Both the source radii and fitted radius 
parameters decrease with increasing tot- On the other 
hand, these radii increase with increasing parton scatter- 
ing cross section as a result of both a larger source size 
and stronger collective expansion |24|| . 
Since Eq.(|h gives 



P±t,0±t, 



(5) 



and i?L n = D x 



lde , the ratio i? ut / -Rsidc contains in- 
formation about the duration of emission and has been 
studied extensively [[| g[. This connection between i? ut 
and the emission duration becomes clearer if we neglect 
the x ou t ~ t correlation term D XovLtt pjj in Eq.(^). In the 
inset of the lower-right panel of Fig. g, we show the ratio 
i? ou t/i? s ido for midrapidity charged pions. It is seen that 
the ratio i? ut/ Rsidc obtained from the emission function 
(solid curve) has a value between 1.0 and 1.3 as in predic- 
tions based on the hydrodynamical model with freezeout 
treated via the hadronic transport model However, 
with the radius parameters extracted from the Gaussian 
fit to the correlation function, the ratio i?out/-Rside be- 
comes much closer to 1, similar to the experimental val- 
ues extracted from the measured correlation function [|| . 
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FIG. 3: x ou t — a; s idc (left panel) and x ou t — t (right panel) 
distributions at freezeout for midrapidity charged pions with 
125 < pt < 225 MeV/c. The curve with open diamonds 
represents (x out ) as a function of t. 

To investigate the reason for the large difference be- 
tween the source radii obtained directly from the emis- 
sion function and from the Gaussian fit to the correlation 
function, we show in Fig. |^ the ir ou t — £ S idc distribution 
(left panel) and the x ou t — t distribution (right panel) 
at freezeout for midrapidity pions with 125 < pt < 225 
MeV/c from the AMPT model with string melting and 
Up = 10 mb. It is seen that the emission source is shifted 
in the direction of the pion transverse momentum, i.e, 
(^out) > 0. This positive shift in x out results from the 
collective expansion of the emission source |l4|] . 

If the width in a; out were much smaller than the average 
value (x ou t), the emission source would have a shell-like 
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shape, similar to the emission source with strong opacity 
in the hydrodynamical model jnj. The emission source 
also shows a large halo around a central core. The halo 
consists not only of pions from decays of long-lived res- 
onances such as the w but also of thermal pions. In cal- 
culating the correlation function, we have included pions 
from the decay of 77 resonances as in experiments. Their 
effects on the radius parameters obtained from the Gaus- 
sian fit to the correlation function are thus included. On 
the other hand, we have excluded these pions in evalu- 
ating the source radii from the emission function due to 
the long lifetime of 77. Since long-lived resonances mainly 
affect the correlation function at small relative momenta 
p5| , they are important in determining the A parameter. 

As to the a; ou t — t distribution of the emission func- 
tion shown in the right panel of Fig. |3[ it has a strong 
positive x ou t — t correlation as clearly seen from the solid 
curve with open diamonds, which shows that the aver- 
age value (Xout) increases with the freezeout time t. The 
a^out — t correlation leads to a large positive value for the 
a^out — t correlation term D Xout ^ ±t in Eq.Q. For pions 
included in generating Fig. the value of D Xout ^ x t is 
168 fm 2 and is appreciable compared to 185 and 431 fm 2 , 
respectively, for the first and last terms in Eq.(||). This 
makes it difficult to extract information about the dura- 
tion of emission from the ratio i?out/-Rsidc- Similar results 
have been seen previously in studies based on the RQMD 
model at SPS |^6|, ^7). We note that the imaging method 
psj , developed for extracting the emission function of a 
source from the correlation function, will be very useful 
for verifying the non-Gaussian features of the emission 
source in high energy heavy ion collisions. 

In summary, we have studied the pion interferometry 
in relativistic heavy ion collisions at RHIC. Using a multi- 
phase transport model that includes both initial partonic 
and final hadronic interactions, we find that the two-pion 
correlation function is sensitive to the parton scattering 
cross section, which controls the density at which the 
parton-to-hadron transition occurs in the AMPT trans- 
port model. To reproduce the measured correlation func- 
tion in central Au+Au collisions at y/s = 130 A GeV 
requires both the melting of initial strings to partons 
and a large parton scattering cross section, i.e., a low 
parton density for the transition from the partonic to 
the hadronic matter. We further find that the emis- 
sion source is non-Gaussian in space and time. It not 
only shifts significantly to the direction along the pion 
transverse momentum but also has a strong correlation 
between this displacement and the freezeout time. Con- 
sequently, the source radii extracted directly from the 
emission function are about a factor of 2 to 3 larger than 
the radius parameters extracted from a Gaussian fit to 
the three-dimensional correlation function. Furthermore, 
the ratio i?out/-Rsidc obtained from the emission func- 
tion is larger than that extracted from a Gaussian fit 
to the correlation function, which is found to be close 



to one. Although the correlation function requires only 
the space-time information of pions at freezeout, it is 
shown to be sensitive to the partonic dynamics during 
the early stage of heavy ion collisions. The study of pion 
interferometry thus helps to confirm the formation of the 
partonic matter at RHIC and to study its properties. 
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